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Abstract — In this paper, we propose an approximate semi- 
definite programming framework for demodulation and equal- 
ization of non-coherent ultra-wide-band communication systems 
with inter-symbol-interference. It is assumed that the commu- 
nication systems follow non-linear second-order Volterra mod- 
els. We formulate the demodulation and equalization problems 
as semi-definite programming problems. We propose an ap- 
proximate algorithm for solving the formulated semi-definite 
programming problems. Compared with the existing non-linear 
equalization approaches, such as in [I], the proposed semi-definite 
programming formulation and approximate solving algorithm 
have low computational complexity and storage requirements. 
We show that the proposed algorithm has satisfactory error 
probability performance by simulation results. The proposed 
non-linear equalization approach can be adopted for a wide 
spectrum of non-coherent ultra-wide-band systems, due to the 

I fact that most non-coherent ultra-wide-band systems with inter- 
symbol-interference follow non-linear second-order Volterra sig- 

' nal models. 

I. Introduction 

Ultra-Wide-Band (UWB) communication systems have at- 
tracted much attention recently. The UWB communication sys- 
tems have many advantages including multi-path diversity, low 
possibilities of intercept and high location estimation accuracy. 
However, UWB systems also present many challenges com- 
I pared with narrow-band communication systems. Especially, 
■ the communication channels are frequency selective with a 
large number of resolvable multi-paths. Accurate estimation 
of channel impulse responses is complex and difficult. 
' Existing modulation schemes for UWB can be roughly 
classified into two categories, coherent modulation schemes 
and non-coherent modulation schemes. The coherent schemes 
include direct-sequence UWB and multi-band UWB [2] [3] 
[4]. In these schemes, the demodulation usually depends on 
accurate estimation of channel impulse responses. The coher- 
ent schemes can achieve higher transmission rates. However, 
their complexity and cost are usually high. 

Unlike the coherent modulation schemes, in the non- 
coherent UWB modulation schemes, the demodulation usu- 
ally does not depend on full knowledge of channel impulse 
responses. Therefore, the difficulty of channel estimation is 
largely avoided. The non-coherent schemes include various 
differential encoding schemes, and energy detection based 
schemes (see for example [5] [6] [7]). 

One difficulty with the non-coherent modulation schemes 
is that the signal models are non-linear, if there exists Inter- 
Symbol-Interference (ISI) in the systems [8] [7]. The existing 



linear equalization approaches generally do not work well for 
such non-linear ISI. Because the UWB channels usually have 
long delay spreads, the approach that increases the spaces 
between symbols to avoid ISI, severally limits the achievable 
rates, and therefore is not realistic. 

In [1], a new non-linear equalization scheme based on Semi- 
Definite Programming (SDP) has been proposed. It is shown 
that even though the SDP relaxation approach is sub-optimal, 
the performance loss is usually negligible. In [1], an off-the- 
shelf general-purpose algorithm is adopted to solve the SDP 
programming problems. 

However, general-purpose SDP solving algorithms may 
not be suitable choices for the UWB demodulation and 
equalization scenarios. First, the general-purpose algorithms 
are usually designed to obtain very accurate optimization 
solutions. While, in the UWB demodulation and equalization 
scenarios, only approximate solutions are needed to ensure low 
demodulation errors, because the SDP optimization solutions 
are only intermediate results. By relaxing the requirement on 
the accuracy of optimization solutions, the computational com- 
plexity can be greatly reduced. Second, the general-purpose 
SDP solving algorithms do not utilize problem structures. In 
fact, the computational complexity can be largely reduced by 
utilizing the structure of the problems. 

In this paper, we propose a new iterative algorithm for solv- 
ing the SDP programming problems. The proposed algorithm 
has low computational complexity and storage requirements, 
which make it an attractive choice for low-complexity high- 
speed implementations. First, the algorithm can achieve a 
close approximate solution of the optimization problem after 
only a few iterations. Second, during each iteration, only one 
optimization problem with much smaller problem size needs 
to be solved. More precisely, the problem size is equal to the 
number of bits in one signal block, while, the problem size of 
the original matrix optimization is proportional to the square 
of the number of bits in one signal block. The correctness 
and convergence of the algorithm is proven in this paper 
We also show by simulation results that the demodulation 
and equalization algorithm has satisfactory error probability 
performance. 

In this paper, we demonstrate the performance of the 
proposed non-linear demodulation and equalization scheme 
on differential UWB systems. In fact, the proposed algorithm 
can also be applied on other non-coherent UWB systems, 
because many non-coherent UWB systems have the same 
non-linear second-order Volterra signal models. One thing we 
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wish to stress is that certain channel parameter estimation is 
needed in the proposed demodulation algorithm. However, the 
estimated model is at the symbol level, rather than at the 
Nyquist frequency level. The complexity of this partial channel 
estimation is acceptable. 

The rest of this paper is organized as follows. In Section HIl 
we describe the signal model. The SDP problem formulation 
is presented in Section |lll] We present the proposed demod- 
ulation and equalization algorithm in Section |IV] Numerical 
results are presented in Section |V] Conclusions are presented 
in Section |VT] 

Notation: we use the symbol <S to denote the set of symmet- 
ric matrices. Matrices are denoted by upper bold face letters 
and column vectors are denoted by lower bold face letters. 
We use A ^ to denote that the matrix A is positive semi- 
definite. We use a > to denote that the elements of the 
vector a are non-negative. We use Ai_j to denote the element 
of the matrix A at the i-th row and j-th column. We use a; 
to denote the i-th element of the vector a. We use A*^ and 
a* to denote the transpose of the matrix A and the vector a 
respectively. We use tr{A) to denote the trace of the matrix 

A. We use A B Xo denote the inner product of matrices A and 

B, that is A ■ B ~ tr{A*-B). The function sign(-) is defined 
as. 



sign(2;) 



1, ifx>0, 
— 1, otherwise. 



(1) 



II. Signal Model 



Information Bearing Signal 



where 'w{t) is the transmitted pulse, ai[n] is the pulse polarity 
for the i-th pulse of the n-th symbol, ti [n] is the pulse time for 
the i-th pulse of the n-th symbol. Each block has symbols, 
and each symbol corresponds to Np pulses. 

Denote the data symbol by d[n] G {—1,+!}. The data 
symbols are differentially encoded as. 



ai\n\ 



aArp_i[n - l]d[n - l]6jVp-i, if i = 
ai^i[n]d[n]bi^i, otherwise 



(4) 



where, foo, ^i, ■ • • j ^Wp-i is the pseudo-random amplitude code 
sequence, hi G { — 1,+!}. The pulse time 



ti [n] = nTs 



(5) 



where Tg is the symbol duration, Ci is the relative pulse timing. 

The pilot signal (t) is introduced to facilitate timing syn- 
chronization and partial channel estimation. Guided intervals 
are introduced between blocks of information bearing signals 
and pilot signals, so that all inter-block-interference is avoided. 

Similarly as in [1], at the receiver side, an auto-correlation 
receiver is used. Denote the received signals corresponding 
to one block of information bearing signals by z[m], m = 
1,2, Nr- The signal model of the system is a second-order 
Volterra model as follows. 

z[m] = (r + PdYQ*B[m]Q{r + Pd) + noise terms, (6) 

where Q, P, r are constant matrices and vectors , and B[n] are 
matrices that depends on the wireless channel (more detailed 
definitions can be found in [1]). We assume that the matrices 
B[m] can be estimated accurately by using the pilot signals. 




Pilot Signal 



Fig. 1. signal is transmitted in a bloclc by block fashion 

In this paper, we consider the differential UWB systems. 
We assume that information is transmitted in a block by block 
fashion as shown in Fig.[T] That is, the transmitted signal s(t) 
can be written as, 

oo oo 

<t) = 'kit ~ m) + ^ slit - m - Tfc) (2) 

fe=0 k=0 

where s\.{t) is the signal waveform for the fcth block of 
information bearing signals, and s^,(t) is the waveform for 
the fcth block of pilot signals. 

The waveform for one block of information bearing signals 
can be written as, 

^k{t)^ £ a^[n]w{t~t,[n]) (3) 

n=0 1=0 



III. SDP Problem Formulation 

Similarly as in [1], we reformulate the difficult discrete 
optimization problem into a matrix optimization and relax it 
into an SDP problem. The SDP formulation in this paper is 
slightly different from the one in [1]. Instead of introducing 
auxiliary variables, we formulate the SDP problem with the 
following convex objective function f(U). 

f{U) = Y {-^H - r*Q*B[m]Qr - r*Q*B[m]QPd 

m—l 

- r*Q^B[mYQPd - tr {DP^Q^ B[m]QP}Y , (7) 

where U is, a. + 1 hy Ni, + 1 positive semi-definite 
symmetric matrix, D denote the sub-matrix of U formed by 
selecting the last rows and columns, and d is a vector 

d= [t7i,2,...,t/i,A,, + i]* 

The convex SDP problem is summarized as follows. 

min fiU) 

subject to: Un,n = 1, for all n, (8) 
UeS, (9) 
U hO. (10) 

The demodulation result is obtained from the solution of the 
above SDP problem by thresholding. That is, the demodulation 
result for the nth symbol is obtained as sign(J7i.„+i). 
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IV. Approximate Semi-definite Programming 
Algorithm 

In this section, we propose a new approximate algorithm 
of solving semi-definite programming. The algorithm is a 
generalization of Hazan's algorithm on approximate semi- 
definite programming [9]. Hazan's algorithm considers a spe- 
cial class of SDP optimization problems, where the constraints 
are total trace constraints. Such SDP optimization problems 
usually arise in Quantum State Tomograph (QST) problems. 
The algorithm proposed in this paper considers the class of 
problems with the constraints that the diagonal elements of 
the matrix must be one. 

We consider the following SDP optimization problem. 

mmf{X) 

subject to: diagonal elements of X are zeros, 
X is symmetric, 

X+IhO, (11) 

where, X is a square matrix, / is the identity matrix with 
the same numbers of rows and columns. Without loss of 
generality, we assume that /(•) is independent of the diagonal 
elements of the matrix X. We also assume that /(•) has a 
bounded curvature constant C / . The curvature constant C / is 
defined as follows. 

Cf = sup ^ [fiY) - fix) + {Y- Xy V fix)] (12) 

where, X +1 h 0,Z + 1 h 0,F X + - X), and 
all diagonal elements of X, Y, Z are zeros. Clearly, the 
convex-SDP optimization problem in the previous section can 
be reduced into the above optimization problem and solved. 

Before going into details of the proposed algorithm, we need 
some basic facts on matrices. These facts will be presented in 
Section IIV-AI The dual function of /(X) will be discussed 
in Section IIV-BI The algorithm will be presented in Section 
IIV-CI The correctness and convergence of the algorithm will be 
proved in Section IIV-DI Certain discussions will be presented 
in Section HV^ 



Lemma 4.2: Let Xi X2 be two symmetric matrices, such 
that the smallest eigenvalues of the matrices are greater than 
-1, 



Amin(-X'i) > —1, A,nin(-X^2) > ^ 1- 



(14) 



Let X be a linear combination of Xi and X2- That is X = 
PXi + (1 - I3)X2, where < /3 < 1. Then, the smallest 
eigenvalue of X is also greater than —1, 



Amin(A:) > -1. 



(15) 



Proof: 



Aminf-X") = min v Xv 
..t 



= min v\pXi + (1 - (3)X2)v 
ll"ll=i 

>/?(-!) + (l-/3)(-l) = -L 



(16) 



B. Weak Duality 

The proposed algorithm is based on iteratively reducing the 
duality gap between the primal function and its dual function. 
For a primal function /(X), we define the dual function wiX) 
as 



wiX) = max wiX,X) 



max fix) 
vf{x)+xyo 



-X- v/(X)-ir(A), (17) 



where, A — diag(Ai, A2, . . . , An) is a diagonal matrix. 

Theorem 4.3: {Weak Duality) Denote the minimizer of the 
optimization problem in Eq. [TT] as X* . Let X be a feasible 
point. Then, the following weak duality inequalities hold. 

/(X) > /(X*) > wiX) (18) 
Proof: Given a function /(X), the corresponding La- 
grangian function can be written as 



/(X)-y-(7 + X)+A-X 



(19) 



where y is a symmetric positive semi-definite matrix. 

We can rewrite the function /(X*) in a min-max form as 
follows. 



A. Some Basic Facts 

Lemma 4.1: Let X be a symmetric matrix with all the 
diagonal elements being zero. Then / + X is positive semi- 
definite, if and only if Aniin(X) > — 1, where Aniin(X) denotes 
the smallest eigenvalue of X. 

Proof: Necessary condition: assume that X is positive 
semi-definite, then 

Amin(X) — mill v*Xv, 

ll«ll=l 

= min (7 + X)v - v*^Iv, 
ll«ll=i 

> min Q~v*Iv = -1. (13) 

Sufficient condition: It is sufficient to show that t)*(7 + 
X)v > for all v with \ \v\ \ = 1. The above statement follows 
from the fact that v^Xv > Aniin(X)||t;|p > -1. ■ 



/(X*) = min /(X) 



= mm 
xes 



max [/(X) - y • (/ + X) 
vya.x 



XX] 



(20) 



This is because 
max [/(X) 

fiX), 



-00, 



V-iI + X)+X-X] 
if X hO and 

diagonal elements of X are zeros 
otherwise 



(21) 



By the max-min inequality (see for example, [10] page 238, 
Eq. 5.47), we can lower bound /(X*) as follows. 



fix*) = min 



max [/(X)-y-(7 + X) 
vyo,\ 



XX] 



> max 

V>-0,A 



min IfiX) 
xes ^ ^ ' 



V^- (7 + X) + A-X] 



(22) 
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Let us assume that Vq and Aq are symmetric and diagonal 
matrix respectively, such that the following equations hold for 
a feasible point Xq. 



V/(J^o)-Vo+Ao=0, 
V/(Xo)+Ao hO. 

By the above discussions, we have 

}{X*) 



(23) 
(24) 



> max 



min [f{X) - V -{I + X)+\-X] 



> min [f{X)-Vo-{I + X)+Xa-X], 

f{Xo)-Vo-{I + Xo)+Xo-Xo, 

^=V(Xo) - (v/(^o) + Ao) ■ (/ + Xo) 
= f{Xo)-vf{Xo)-X -Xo-I, 



^0 ■ Xq, 



(25) 



where, (a) follows from the fact that Xq is exactly the mini- 
mizer, and (b) follows from the definition of Vq. Therefore, 



fix*) > max f{Xo) - Xo ■ V/(^o) - HX), 



> w{Xo). 



(26) 



The theorem follows from the fact that Xq, Vq, and Ao are 
arbitrary. ■ 
The above weak duality theorem provides a way to estimate 
how far a feasible point X is away from the optimal solution. 
We define 



h{X) = f{X) - f{X*) 
g{X) = f{X) - w{X) 



(27) 
(28) 



By the weak duality theorem, we have h{X) < g{X). 

In order to evaluate the dual function w{X), the following 
optimization problem needs to be solved. 



min Xi 

i 

subject to 7i > 0, for all i 



(29) 



where 7^ is the ith eigenvalue of the matrix \/ f{X) + A. 

Lemma 4.4: Let 7^ denote the ith eigenvalue of the matrix 
\/f{X) + X. Let Vi denote the corresponding eigenvectors. 
Then, 



(30) 



where, Aji and AA^ are the infinitesimal differences, Vij 
denotes the jth element of the vector Vi. 

Proof: It is clear that there exists a decomposition of 

vfix)+x. 



v/W + A = yAF* 



(31) 



such that y is a unitary matrix and A is a diagonal matrix. 
In fact, V ~ [vi,V2, ■ ■ ■ ,Vn] and A — diag(7i) is a such 
decomposition. 



Let AV and AA be the corresponding infinitesimal differ- 
ences of V and A respectively. Then, we have 



{V + AV){A + AA){V* + AV*) = V/(X) + A 



AA, 
(32) 



{V* + AV*){V + AV)=I. (33) 

From Eq. [33] and the fact that V is unitary, we have 

V*AV + AV*V = 0. (34) 

Since V*AV and AV*V are the transpose of each other, 
we conclude that the matrices V*AV and AV*V are anti- 
symmetric and their diagonal elements are all zeros. 
From Eq. |32] we have 



VAAV* + AVAV* + VAAV* AA. 



(35) 



Multiplying the above equation by the matrix V* at the left 
side and the matrix V at the right side, we obtain 



AA + V*AVA + AAV*V = V*AXV. 



(36) 



Since the diagonal elements of the matrices V*AV and AV*V 
are all zeros, the diagonal elements of the matrices V*AVA 
and AAV*V are also zeros. Therefore, we conclude that the 
diagonal elements of AA and V*AXV are identical. The 
theorem then follows from the fact that the ith diagonal 
element of the matrix F*AAV is J2ji''^ij)^AXj. ■ 
Lemma 4.5: In the optimization problem in Eq. |29l Let A* 
be the minimizer Let 7^ denote the ith eigenvalue of the matrix 
\/f{X) + A*. Let Vi denote the corresponding eigenvectors. 
Then, there exist a set T C {1, 2, . . . , n} and a vector y, such 
that 

y>0, (37) 
V*y=[l,...,l]\ (38) 

vl (v/(X) + A*)t;, = 0, for all ^ G T, (39) 

where n is the number of rows of matrix A, is a matrix 
such that each row of V is [vfA for one i ^ T. That is. 



V 



(40) 



ifc,2 ■ ■ ■ "ik,n 

where ii,...,ik £ T. 

Proof: Due to the nature of the optimization problem, 
there exist at least one active constraint at the minimizer. 
We say that an inequality constraint is active at a feasible 
point, if the inequality constraint holds with equality. In this 
optimization problem, the ith inequality constraint is active, 
if 7i = 0. Let T denote the set of indexes of all active 
constraints. Then, for all i eT, 7; = 0, 



v\ (v/(X)+A*)u, = 0. 



(41) 



Due to the Karush-Kuhn-Tucker (KKT) Theorem ( see [11] 
Theorem 20.1 . Page 398), there exists a vector y such that 



2/>0, 



(42) 
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Y^y^Vl^^ V^I = [1, 1, • • ■ , 1] 



(43) 



According to Lemma 14.41 V7j = bfii '^'i2T ■ ■ ^ '^'inY ■ There- 
fore 



(44) 



The lemma follows. ■ 
Corollary 4.6: For i <E T, define ai — v^X/ f{X)vi. Then, 

E^»* = -E2^«"' ^45) 

Proof: 

^ A* = [1, . . . , l]diag(A*) = j/V(diag(A*)) 

i 

t 



= y 



[...,Vi'X*v^,...y 



(46) 



where, diag(A*) denote the column vector that consists of 
diagonal elements of A*, and (a) follows from Eq. [39l ■ 

C. The Algorithm 

The proposed algorithm is summarized as follows. 

• Step 1: set k=l, set X to a feasible point; 

• Step 2: calculate the gradient \/f(X); 

• Step 3: solve the optimization problem in Eq. |29l obtain 

yi, Vi for i e T; 

• Step 4: calculate the function g{X), if g{X) is less than 
a certain threshold, go to step 8, otherwise, go to the next 
step; 

• Step 5: update AX as follows. 



X I 



where, (3k is a predefined step size parameter; 
. Step 6: update X = X + AX; 

• Step 7: set k=k+l, go to step 2; 

• Step 8: return X, stop. 

D. Correction and Convergence 

In this subsection, we show that the proposed algorithm is 
correct and converges. 

Theorem 4.7: In the proposed algorithm, the diagonal ele- 
ments of X are zeros, and Aniin(-X^) > ^ 1- 

Proof: Prove by induction. It is sufficient to show that 
X + AX satisfies the above conditions, if X satisfies the 
conditions. 

Note that the fcth diagonal element of the matrix 

J^iGTVi'^M is equal to Y^it^r Vi'^'ik' ^l^o equal to the 
fcth element of V*y. From Lemma 14.51 we have that the fcth 



element of V^y is one. Therefore, the diagonal elements of 
the matrix J^ieT Vi'^i'i^i — X — 7 are all zeros. The diagonal 
elements of the matrix X + AX are also all zeros. 

We can show that A,nin(X + AX) > — 1, if we can show 
that 



\ieT J 



A,„i„ (X) > -1. 



> -1, 



(48) 



(49) 



This is because of Lemma 14.21 and X + AX being a linear 
combination of the above two matrices 



X + AX - /3fc 



\ieT 



ViViVi 



{l-(ik)X (50) 



Eq.|49l follows from the given hypothesis. Eq. [48] follows from 
Lemma 14.11 and X^ier yi'^i''^i being positive semi-definite. 
The theorem is proven. ■ 
Theorem 4.8: In the proposed optimization algorithm, let 
Xfe denote the value of X after k iterations. Then, the gap 



hiXk+i) < (1 - Pk)h{Xk) + PlCf. 



(51) 



Therefore, h{Xk) goes to zero, and /(X^) goes to /(X*), 
for properly chosen step size parameters f3k- 

Proof: First, we wish to show that the following equality 
holds. 



(52) 



vier 



(47) rpjjg reasoning is as follows. 



= Y.{tr {y.,v,vlx7.f{X))) 

'^Y.y.{tr{vlvfiX)v.)) 
= Y.y^{vlvf{X)v..) 



(fc) 



E yi^i 



(53) 



where, (a) follows from the the property of trace, tr{AB) — 
tr{BA) for all matrices A and B, and (b) follows from the 
definition of ai. 
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The value of f{Xk+i) can be upper bounded as follows. 

< fiXk) + f3k yrvM - Xk - ■ vfiXk) + PlCf 




f{Xk) + Pk v^o., - Xk ■ VfiXk)^ + f3lCf 



= f{Xk)-f3k9iX,) + PlCf, 



(54) 



where, (a) follows from the definition of C/, (b) follows from 
the fact that the diagonal elements of sjf{X) are all zeros, 
/ • Sjf{X) — 0, (c) follows from Eq.|52l and (d) follows from 

Eq.m 

Therefore, we have 

h{Xk+i) = f{Xk+i)-.nx*) 

< /(Xfc) - f{X*) - hg{Xk) + PlCf 
<h{Xk)-Pk9{Xu)+PlCf 
<h{Xk)-PkhiXk) + l3lCf 
<{l-l3k)h{Xk)+l3lCf. (55) 

The theorem follows. ■ 
E. Discussion 

One character of the proposed algorithm is that a close 
approximate solution can be found after only a few iterations. 
By Theorem 14.81 we can see that the optimal step size 
parameter /3fe at the kt\\ iteration depends on the current gap 
h{Xk) and C/. At the first several iterations, the parameter Pk 
can take larger values, and the gap h{Xk) decreases quickly. 

Because the solution of the SDP optimization problem is an 
intermediate result in the demodulation and equalization algo- 
rithm, approximate solutions are usually sufficient to ensure 
that the demodulation results are correct with high probability. 
In fact, we find that only few iterations are usually needed 
to ensure low demodulation error probability by simulation 
results. 

During each iteration, one optimization problem needs to 
be solved to calculate the dual function. However, compared 
with the original matrix optimization problem with approxi- 
mately optimization variables, the optimization problem in 
dual function calculation only has n optimization variables. 
Therefore, the optimization problem in each iteration can 
be solved with lower computational complexity and storage 
requirements. 

Overall, the proposed algorithm has lower computational 
complexity and storage requirements. It is an attractive choice 
for high-speed real-time demodulation implementations. 



V. Numerical Results 

In this section, we present simulation results for the pro- 
posed demodulation and equalization scheme with approx- 
imate SDP programming. We assume that the transmitted 
pulses are the second derivative Gaussian monocycles. 



w{t)= l-47r(t/r„)^ exp I -27r (t/T„i)^ I 



(56) 



where r„j — 0.2877 nanosecond. Each information bearing 
signal block consists of Nb = 10 symbols and each symbol 
corresponds to Np = 4 pulses. The symbol duration T, = 8 
nanoseconds. 

We use the IEEE 802.15.4a channel models as described in 
[12]. Two types of channel models CMl and CM6 are used 
for simulation. We illustrate the bit error probability of the 
proposed demodulation and equalization algorithm in the case 
of CM6 channel models in Fig. |2] A typical channel impulse 
response of the CM6 model is shown in Fig. [3] We illustrate 
the bit error probability of the proposed demodulation and 
equalization algorithm in the case of CMl channel models 
in Fig. m A typical channel impulse response of the CMl 
model is shown in Fig. |5] The numerical results show that 
the proposed demodulation algorithm has satisfactory bit error 
probabihty performance. 



Bit error probability 




4 5 6 7 8 9 10 11 12 13 
energy per bit to noise power spectral density ratio in db 



Fig. 2. Bit eiTor probabilities for the CM6 channel model. The X-axis shows 
energy per bit to noise power spectral density ratio Ei,/No in dB 



VI. Conclusion 

In this paper, we propose an approximate semi-definite 
programming framework for demodulation and equalization 
of non-coherent UWB systems with inter-symbol-interference. 
The proposed algorithm has low computational complexity and 
storage requirements, which make it an attractive choice for 
real-time high-speed implementations. Numerical results show 
that the proposed approach has satisfactory error probability 
performance. The proposed approach can be adopted in a wide 
spectrum of non-coherent UWB modulation schemes. 
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Channel impulse response 
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Fig. 3. Typical cliannel impulse response in the CM6 channel model 



Bit error probability 




8 10 12 14 16 18 

energy per bit to noise power spectral density ratio in db 



Fig. 4. Bit error probabilities for the CMl channel model. The X-axis shows 
energy per bit to noise power spectral density ratio Ei,/No in dB 
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Fig. 5. Typical channel impulse response in the CM6 channel model 
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